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A theory is developed that gives the diffusion coefficient in strained systems as an exponential function of 
the strain. This theory starts with the statistical theory of the atomic jump frequency as developed by 
Vineyard. The parameter determining the effect of strain on diffusion is related to the changes in the inter- 
atomic forces with strain. Comparison of the theory with published experimental results for the effect of 
pressure on diffusion shows that the experiments agree with the form of the theoretical equation in all cases 


within experimental error. 

I. INTRODUCTION 

S INCE the diffusion rate in a crystal depends on the 
atomic interaction energy, and since this energy 
depends on the interatomic distances, it is to be ex- 
pected that the diffusion coefficient of a migrating 
species will be altered by a strain superimposed on the 
crystal. Experimental evidence shows that the change 
in the diffusion coefficients resulting from strains can be 
considerable. Uniaxial elastic strain can increase the 
self-diffusion coefficient by as much as a factor of two 1 * * 4 * 
and large hydrostatic pressures may decrease the self- 
diffusion coefficient by as much as an order of 
magnitude . 2-6 

The theory of the effect of pressure on diffusion has 
been examined on the basis of the dynamic theory of 
diffusion . 6 ' 7 In this theory, the pressure effect is repre- 
sented by a parameter that is a function of the normal 
mode vibrations of the atoms in the crystal, and the 
diffusion coefficient is an exponential function of the 
pressure. 

1 T. Liu and H. G. Drickamer, J. Chem. Phys. 22, 312 (1954). 

* Norman H. Nachtrieb, Wright Air Development Center 
Technical Report No. 55-68 (unpublished). 

• J. Petit and N. H. Nachtrieb, J. Chem. Phys. 24, 1027 (1956). 

4 W. Jost and G. Nehlep, Z. physik. Chem. 34, 348 .(1936). 

6 Norman H. Nachtrieb, Henry A. Resing, and Stuart A. Rice, 
J. Chem. Phys. 31, 135 (1959). 

4 Stuart A. Rice, Phys. Rev. 112, 804 (1958). 

7 Stuart A. Rice and Norman H. Nachtrieb, J. Chem. Phys. 31, 

139 (1959). 


The dynamic theory of diffusion was developed as an 
alternative to the absolute rate theory of diffusion, 
since it was believed that the absolute rate theory de- 
pended on the postulate that the jumping atom spends 
a long time at the top of the potential barrier. However, 
it can be shown that the theory of the jump frequency 
can be developed without reference to such a postulate 8 
by considering the motion of a representative point in 
phase space. The jump frequency then depends on the 
rate at which phase points move over the potential 
maximum in configuration space, and not on the length 
of time the phase points spend at the maximum. In 
view of this situation, it is of interest to investigate the 
effect of strain on diffusion in terms of the statistical 
rate theory. 

The statistical rate theory of diffusion in strained 
crystals as developed in this paper shows that the dif- 
fusion coefficient is an exponential function of strain, 
and that the strain effect can be represented by a pa- 
rameter that is a function of the interatomic forces. 
The rate theory, therefore, has an advantage over the 
dynamic theory in two respects: First, the effect of 
strain on diffusion in different materials can be corre- 
lated with the interatomic potential energy, and second, 
the interatomic forces provide a basis on which to 
calculate the magnitude of the strain effect for different 
diffusion mechanisms. Accordingly, the possibility pre- 

4 George H. Vineyard, J. Phys. Chem. Solids 3, 121 (1957). 



sents itself of deciding among alternative diffusion 
mechanisms from a comparison of the results of experi- 
ments on the effect of diffusion in strained systems with 
theoretical calculations. Such a program would be con- 
siderably more difficult in the framework of the dynamic 
theory. 

The general equation for the diffusion coefficient for 
the flow of a single species in an isotropic solid may be 
written 

D=a\ 2 nT, (1) 


JatticesHTm^ffie ^gpfitial energy in Eq. (2) is given by 

<P= <p(qi,taf>), (3) 

where <?,■ represents the set of all atomic coordinates 
and e a p represents the set of six independent strain 
components. 

The potential <p can be expanded as a Taylor series 
in the strains about the point of zero strain with the 
result that 


where D is the diffusion coefficient, X is the lattice 
parameter, n is the concentration of carrier defects, F 
is the jump frequency, and a is a constant that is deter- 
mined by the crystal structure. -In the following sec- 
tions expressions are derived for the effect of homo- 
geneous static strains on the jump frequency and va- 
cancy concentration. The resulting equations are put 
into a form in which comparisons can be made with 
existing experimental data. 

II. DEPENDENCE OF JUMP FREQUENCY 
ON STRAIN 

According to the statistical theory of rate processes, 
the jump frequency is determined by the ratio of two 
configurational integrals, one referring to the activated 
state and the other referring to the normal state. In 
analyzing the effect of strain on the jump frequency, the 
formulation of the rate process theory in solids given by 
Vineyard 8 is used, in which the jump frequency is given 
in terms of these integrals by 


a.e 

^ (4) 

a t (3,p,<r 

where the coefficients C a g and C a $ p , are defined by 


/ d<p\ 

(5) 

Ca0=l I , 

\de a g/ Qj.o 


1/ dV ^ 


appv I ) 

(6) 

2 3j'(0 



The subscripts indicate that the derivatives are evalu- 
ated when the strains are zero and the coordinates 
have the value q,. 

Substituting Eq. (4) into Eq. (2) gives the jump 
frequency in terms of the strain r 

r(e^) 


je-vWdo j J e-^ kT dA, (2) 

where k is Boltzmann’s constant, T is the temperature, 
and <p is the potential energy of the system as a function 
of all the coordinates of all the atoms in the crystal. 
The integral in the numerator of Eq. (2) is evaluated 
over a hypersurface u in the configuration space such 
that the surface passes through the point corresponding 
to the diffusing atom at its activated position with all 
other atoms at their equilibrium positions. The hyper- 
surface is also required to be perpendicular to contours 
of constant potential energy in the configuration space. 
The hypersurface defined in this manner divides the 
configuration space into two symmetric parts. The 
integral in the denominator is evaluated over the config- 
uration volume A of one of these symmetric parts. 

Equation (2) was derived for the case of an unstrained 
crystal. However, it is applicable to strained crystals if 
the potential energy is taken to be a function of the 
six strain components t a g as well as the atomic co- 
ordinates qi. A similar procedure has been used by Born 9 
in an analysis of the statistical mechanics of crystal 



(>—<p(.Qj,0)/kT 



1 

kT a,0 



r e - f («.o)/*r e xp^ ~L £ 0.0 W (7) 
"A V kT a,0 / 


where terms of order higher than the first have been 
ignored. It will_be shown later that the first-order con- 
tribution of the strain to the jump frequency depends 
on the difference of the average value of C a g evaluated 
near the normal configuration and near the acti- 
vated configuration, and on similar differences 
in the averages of C a g pa , etc. It is extremely difficult 
to give an a priori estimate of the relative magnitudes 
of these differences. At any rate, for small enough 
strains the first-order terms predominate and the higher 
order terms can be neglected. It will be seen later that 
the form of experimental results is adequately de- 
scribed by considering only the first-order terms in the 
strains. For zero strain, Eq. (7) gives the jump fre- 
quency as 


r 0 = ai a g- r (.<,i.O)lkT dA (g) 


9 Max Born, Proc. Cambridge Phil. Soc. 36, 160 (1940). 
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Now take the ratio of Eqs. (7) and (8). The result is 

r («->_/ / 1 vr \\ / 

\ expl X CapPaP I / / 

r 0 \ V kTa.p // J 



1 

kT a.P 



The angular brackets indicate that a statistical average 
has been taken of the quantity within the brackets, and 
the subscripts <r and A indicate that the averages are 
taken over the regions of configuration space a and A , 
respectively. The explicit expressions for these averages 
are 


^exp^ X 

\ V kTa.p //, 

= j 6 ^ 0)lhT eXP (“ ~ £ CaPtaP^dc j 

f f-itvMMdv, (10) 

^ ” O 

\ eX Pl H C a ptap I / 

> V kT a,P / / A 

= J* e~ vi "’ fi)lkT exp^— — X C a pi a p^dA j 

f e -*(«i.o)!kT dA ( lla ) 

J A 


Taking logarithms of Eq. (11) and utilizing the fact 
that ln(l — x)~ — x for small x, gives 

T(e a p) 1 

In = — {X (C a p)J]tap), (14) 

To kT a.p 

or, defining a parameter tn a p by 

m a p (C a p) A (Ca^)ff* (lo) 

Equation (14) can be written as 

r(e a/ j) = r 0 expf — X m a pt a p\ (16) 

\kT a.B / 

Since r 0 can always be written as 8 

r = v * e -PEVkT > (17) 

where A E* is the energy of activation for the atomic 
jump and v* is an effective frequency, it is evident from 
Eq. (16) that the strain aSects the jump frequency by 
an effective change in the energy of activation. 

Equation (16) shows that the jump frequency has a 
simple exponential dependence on the strains and that 
this dependence is controlled by the derivatives of the 
potential energy with respect to the strains evaluated 
at the saddle point of the activated state. 

Equation (16) gives the general relation between the 
jump frequency and the strain that will be used in this 
paper. 

To illustrate the application of Eq. (16), three special 
cases will be considered : 

(1) Uniform compression or expansion, in which 

Pxx~ Pyy~ Pzz Pj ( 1 $) 


For small strains and high temperatures, the conditions 
under which the experimental effects of strain on dif- 
fusion are usually determined, the exponents in Eq. (9) 
can be expanded into a series, and only the first two 
terms need be retained. Thus, Eq. (9) can be written as 


T(u t ) 


= ^1 — — X ( C a p)ctc,p ^ j 


( 1 £ ( C„p) A tap\ (lib) 

V kT a.p / 

where ( C a p) c and (C a p) A are given by 

(C a p),= f (—) / 

*9 \ de a p f ff/,o • 

jje-vUi.WkTfo' (12) 

(C a p) A = f (—) e~*W**dA / 

^ A \dtaP' tj.O ' 


X 


e-p(o,MkT dA _ ( 13) 


(all other strains=0) 

(2) Simple shear, in which 

Pxy Pyx Pay ( 19 ) 

(all other strains =0) 

(3) Simple elastic tension or compression in the * 
direction, in which 


11 tL ’ ( 20 ) 

where n is Poisson’s ratio. For these three cases, Eq. 
(16) gives the following results: For uniform compres- 
sion or expansion, 

r(e) = roe (3 ’"'* r >‘, (21) 

where 

m=(C) A -(C)„ (22) 

and 


(cy 


(C),= J (-pj <r*W" kT dG / J e-Wo.o)/*^, (23) 

= /(j") c*' imvr dA j 

J e-*i*,.o)/kT dA ( 24 ) 
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For simple shear, 

r(«.) = r 0 e®-'* r, ‘*, ( 25 ) 

where 

(26) 

( C,)a and (C,)„ are given by statistical averages similar 
to Eqs. (23) and (24), i.e., 


<C.).= 



<&>,-((-) \. 


(28) 


For simple elastic tension or compression in 
direction, 

r(£ Z .)=roe ( ”“ /ir)(I - 2 '‘>, 

where 

mi={Ci) A —(Ci)„ 


or 


”*'=<(-) >-<(-) >• 

o' a ' Vd£L'9>. of » 


the x 


(29) 

(30) 


in. DEPENDENCE OF VACANCY CONCENTRATION 
ON STRAIN 

The n appearing in Eq. (1) has a different significance 
for different mechanisms of diffusion. Broadly speaking, 
n is the probability that a diffusing particle has a site 
available to jump into. For dilute interstitial diffusion 
this probability is nearly unity whether or not the 
system is strained. For diffusion by a vacancy mecha- 
nism, however, n is the atomic fraction of vacancies in 
the crystal, given by 

n=n v /N T , (31) 


vacancy and \po is the energy of the perfect crystal. 
Performing the integrations over the momenta py con- 
verts Eq. (32) into 


n v — J ■ • • Je *" ( q ’ ),kT II dq; XI ( v 0 )k j 

f f e ~ v ° {gi),hT n dq,- n (v v )k, ( 33 ) 


where <p„(q) and <po(q) are the potential energies in a 
crystal containing a vacancy and in a perfect crystal, 
respectively, each taken as a function of all the co- 
ordinates; (v,)k is the frequency of the kth. vibrational 
mode in the crystal containing a vacancy; and (vo)k is 
the frequency of the &th normal mode in a perfect 
crystal. In a strained crystal, the <p and the v must be 
written as functions of strain, so that (33) becomes 

n v (e a fi) 

■■■ f exp[- <p v (qj,e a fi)/kT] II dq/JJ M* 

J j k 

If exp[- <Pa(q 1 ,t«p)/kT~] II dqj I [ (i//)*. 

(34) 

An estimate of the effect of strain on the frequency 
ratios can be made from Griineisen’s relation 10 

d lav/d lnF= — 7 , (35) 

where V is the volume and 7 is a positive constant. 
Integrating Eq. (35) for each vibrational mode as the 
crystal goes from the strained to the unstrained state 




where n v is the vacancy concentration, and Nt is the 
total number of lattice sites per cubic centimeter. It 
is therefore necessary to investigate the variation of n v 
with strain. 

The atomic fraction of vacancies in a crystal at equi- 
librium is given by (see Appendix) 



•=/'/- »(j I f dpydqy j 

f " f e ~ Mpi ' vi>/,cT II dpydqy, (32) 


where 57 is the volume change arising from the strain. 
Griineisen’s relation, therefore, leads to an equality of 
frequency ratios in the strained and unstrained systems : 

n Mk/n M k =u (vo)k / n w*. ( 38 ) 

k k k k 


where \p v is the energy of the crystal containing a Therefore, the ratio of Eqs. (34) and (33) is 

f e5L T’[~ II Mi I J' ■ - • 

xf f e ~V0(qj)/kT n V/ "/ exp[- <pa(qi,(cii)/kT2 II dqy. (39) 


10 John Clarke Slater, Inlroductionto Chemical Physics (McGraw-Hill Book Company, Inc., New York, 1939), p. 238. 
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The potential-energy functions <p v (qi,e a p) and <po(q,e a p) 
can be expanded as Taylor series in the strain just as 
in the development beginning with Eq. (4) and leading 
to Eq. (16). The result is 


where 




= expf — Y it> a ^e a p\ 

\kTa.f> / 


»*«=<(-) >-((-) >■ 

''■dtaP'tj, O' ' ' dt a p/ aj.O' 


d<p v \ . 


( 40 ) 

(41) 


For the special case of uniform compression or ex- 
pansion, Eqs. (40) and (41) become 

n c (e) 

=e (3 „ltTU y ( 42 ) 

n v ° 

-<©)„><©„> » 

where e is the volume strain. 

IV. PROOF THAT THE STATISTICAL AVERAGES 
((d<p/dz)qj, o) DO NOT VANISH 

The preceding theory depends on the statistical aver- 
ages of the derivatives of the potential energy of the 
crystal with respect to strain. It has been assumed that 
these averages are not zero, and that a first-order ex- 
pansion in the strains is therefore adequate for small 
strains. This assumption can be justified by expanding 
the crystal energy in normal coordinates. The sta- 
tistical averages of interest all have the form 

<©>/■ /©r-W 

f fe-^Udqj, (44) 


( ©.)■ ■ ■ - 3 f- / (? “K ■ -ii? ? 


Now perform a coordinate transformation according to 
the following definition : 

(52) 

Then, after a few simple algebraic manipulations, Eq. 
(51) becomes 



£ 


exp(—Uj 2 /2kT)dUj (53) 


where <p is the total potential energy of the crystal as 
a function of all the coordinates q y, and the subscript 
zero means that the derivative is evaluated at zero 
strain. If the <?/ are taken to be the normal coordinates, 
<p can be written to the second order as 

y>= p( 0 )+§ Yi ( 45 ) 


where <p(0) is the potential energy when all the atoms 
are at their mean positions, and the Uj are the normal 
mode frequencies. Differentiating Eq. (45) with respect 
to strain gives 


dtp d<p( 0) 
de de 


fE 


do3j 

at 


( 46 ) 


and, since at zero strain the first term on the right is 
zero, 



do>; 

= Y u,——q, 2 - 
i de 


(47) 


For the purposes of this discussion, e will be taken to 
be the strain corresponding to uniform compression or 
expansion, so that for small strains the volume is given 
by 

V— Fo(l+3e), (48) 


Vo being the volume at zero strain. Introducing the 
Griineisen parameter y, by the relation 


d lnco,-/ d lnF= - 7; , (49) 

where the yj are a set of positive constants, and using 
Eq. (48), Eq. (47) becomes 

(-) =-3ZW?A (50) 

V de / o > 

Substituting Eqs. (45) and (50) into Eq. (44) gives 


) ) II dq* j 

(■■ f exp ( Y u,v) n dq,-. (51) 

J J \ 2 kT i / 


and performing the integrations gives 

(0)- 3kT1 r y ‘- (54) 

Equation (54) shows that the averages of the first 
derivatives are never zero and that these averages are 
proportional to the temperature. 

It is extremely difficult to make any a priori decisions 
concerning the signs of m a $ and w a $ defined by Eqs. 
(15) and (41). Such a decision requires a detailed in- 
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vestigation of the variation of localized normal mode 
vibrations with strain in the vicinity of a defect. How- 
ever, on the basis of general physical considerations, it 
is to be expected that both m al 3 and w afl are positive. 

V. EFFECT OF PRESSURE ON DIFFUSION 
CONSTANT 

Using Eqs. (1), (21), and (42), and the fact that the 
lattice parameter in the strained system is (1+e) times 
the lattice parameter in the unstrained system, the 
relation between the diffusion coefficients in the strained 
and unstrained systems for uniform compression or 
expansion is 

Z?( t ) = H u (l+e) 2 exp^— 3ej, (55) 

where D u is the diffusion coefficient in the unstrained 
system, and M is given by 

M=m+w (vacancy mechanism), (56) 

M=m (interstitial, ring or 

exchange mechanism). (57) 

In terms of the volume strain, t=\(AV/Va) for small 
strains, where AF is the initial volume, so that Eq. (55) 
takes the form 

/AF\ / AF\* [M /AF\1 

Kf) =d -( 1 + f) -'bUJ (58) 

Therefore, it is evident that a plot of ln(~D(AF/Fo) 
X (1+AF/Fo) - *] against AF/F 0 should be linear with 
a slope a given by 

a=M/kT, (59) 

and an intercept given by lnU u . 

Several investigators have obtained data on the 
variation of the diffusion coefficients with pressure that 
is suitable for testing Eq. (58). Reference 2 presents 
data for the self-diffusion coefficient as a function of 
pressure for sodium, phosphorous, and mercury up to 
pressures of 12 000, 4000, and 8000 atmospheres, re- 
spectively. The self -diffusion coefficient of liquid gallium 
up to pressures of 10 000 atmospheres is given in 
reference 3. The self-diffusion coefficient for single- 
crystal zinc up to pressures of 10 000 atmospheres for 
diffusion in the directions parallel to and perpendicular 
to the c axis is given in reference 1. The self-diffusion 
coefficient of lead up to pressures of 8000 atmospheres 
at two temperatures is given in reference 5. 

The electrical conductivities of silver chloride and 
silver bromide have been measured as a function of 
pressure up to 300 atmospheres. 4 Since in these silver 
halides it has been shown that the conductivity de- 
pends almost entirely on the mobility of the silver ion, 
the conductivity is proportional to the diffusion co- 
efficient of the silver ion by the Nernst-Einstein rela- 


tion, and the data of reference 4 are suitable for testing 

Eq. (58). 

Plots of the variation of the quantity of log[Z>(AF/Fo) 
X (l+AF/Fo) -i J against AF/F 0 for the self-diffusion 
of sodium, phosphorous, mercury, gallium, and lead 
are shown in Fig. 1. The quantities log[7D (AF/Fo) 
X (l+AF/Fo) -} ] for single-crystal zinc were plotted 
against the fractional change in lattice parameter AX/A, 
since this is a more natural unit for discussing diffusion 
in anisotropic crystals and the linear compressions 
perpendicular and parallel to the c-axis are available. 
The zinc data are plotted in Fig. 2. 

Figure 3 gives log(l/R) plotted against AF/F 0 for 
silver chloride and silver bromide, where R is the re- 
sistivity. The volume change AF/Fo is small enough 
for the pressure range considered so that (1+AF/Fo) - * 
does not appreciably affect the results and can be 
ignored. 

Compressibility data 11-16 were used to obtain the 
appropriate value of AF/F 0 for zinc, sodium, mercury, 
lead, silver chloride, and silver bromide. For gallium, 
AF/Fo was computed from the data of Richards and 
Boyer 17 assuming that the form of AF/Fo as a function 
of pressure is the same as that for mercury. The values 
of AF/Fo for white phosphorus were computed from 
data in reference 16 assuming that the variation of the 
fractional volume change with pressure has the same 
form as that observed 18 for black and red phosphorus. 

In all cases, the available compressibility data were 
extrapolated to the diffusion temperature. 

The linearity of the plots presented in Figs. 1 to 3 
shows that the form of Eq. (58) is valid for those sys- 
tems investigated within the probable inaccuracies of 
the experiments and the calculations. 

The slopes of the plots are given in Table I, where 
a=M/kT and aT are shown for the various materials. 

The fact that aT is so much smaller for the liquid 
metals than for any of the solids including sodium is 
indicative of the difference in the mechanism of dif- 
fusion in liquids and solids. In a liquid, the atoms are 
not constrained to remain at lattice positions, so that 
diffusion occurs by a cooperative process involving the 
migrating atom and its nearest neighbors. Thus, the 
change in the interatomic forces can be kept to a mini- 
mum throughout the diffusion process, and consequently 
aT would be very low. 

From Eqs'. (54) and (59) it is seen that a should be 
temperature independent. For the self-diffusion of lead 
for which pressure data are available at two tempera- 
tures, the value of a is reasonably constant. 

11 P. W. Bridgman, Proc. Am. Acad. Arts Sci. 60, 305 (1925). 

. 12 P. W. Bridgman, Proc. Am. Acad. Arts Sci. 76, 71 (1948). 

13 P. W. Bridgman, Proc. Am. Acad. Arts Sci. 47, 347 (1911). 

14 P. W. Bridgman, Proc. Am. Acad. Arts. Sci. 58-59, 166 (1924). 

15 P. W. Bridgman, Proc. Am. Acad. Arts. Sci. 74, 21 (1940). 

18 P. W. Bridgman, Proc. Am. Acad. Arts. Sci. 62, 207 (1927). 

17 Theodore W. Richards and Sylvester Boyer, J. Am. Chem. 
Soc. 43, 274 (1921). 

18 T. W. Richards, J. Am. Chem. Soc. 37, 1643 (1915). 
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Fig. 1. Variation of log£ZL(l-f A V/V 0 ) - *] plotted against 
volume change (AV/K 0 ) for self-diffusion of various elements, 
(a) Sodium at 363 °K [see enclosed graph], (b) White phosphorous 
at 314°K. (c) Liquid mercury at 303°K. (d) Liquid gallium at 
303°K. (e) Lead at 526.2°K— o ; Lead at 574.2°K— □. 
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(a) 



(b) 


Fig. 2. Variation of logCD.O+AF/Fo) - *] plotted against frac- 
tional change in lattice parameter AX/X 0 for self-diffusion in zinc, 
(a) Zinc at 580°K, perpendicular to c axis, (b) Zinc at 580 °K, 
parallel to c axis. 




Fig. 3. Variation of log(l/2?) plotted against volume change 
AF/Fo for mobility of silver at 573°K. (a) Silver chloride, (b) Sil- 
ver bromide. 


Activation Volume 

The activation volume is ordinarily calculated from 
the relation 


AF* = 


I'd (AG) 1 fd[ln(.D/aXV)]j 


nri, L 


dP 


J T 


(60) 
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Table I. Comparison of values of aT for various solids. 


Solid 

Temperature, 

T 

(°K) 

a 

aT 

Sodium 

363 

27.9 

10 130 

Zinc (X) 

580 

86.6 

50 230 

Zinc (II) 

580 

34.7 

20 130 

Mercury (liquid) 

303 

6.5 

1970 

Gallium (liquid) 

303 

6.5 

1970 

Silver in silver chloride 

573 

91.9 

52 660 

Silver in silver bromide 

573 

128.0 

73 340 

Lead 

526 

139.9 

73 590 

Lead 

574 

121.5 

69 740 


where P is the pressure and A G refers to the Gibbs 
free-energy changes for vacancy formation and for the 
formation of the activated state configuration. This 
free-energy change is calculated from the measured 
diffusion coefficient as a function of pressure. It follows 
from Eq. (58) that (60) may be written 


AF*=-«|- 


[(M/kT)d(AV/V 0 ) 


dP 


■ T 


Since 



d(AF/Fo) 

dP 


(61) 


d(AF/F 0 ) 

dP 



where /3 is the compressibility, the activation volume 
defined by (61) can be calculated from the simple 
formula 


A V*=apkT. (62) 


Table II presents values of the activation volume 
calculated from Eq. (62) at atmospheric pressure for 
those systems for which data are available. 


CONCLUSIONS 

A statistical mechanical theory was developed that 
relates the diffusion coefficient to strain in terms of the 


atomic properties of the system. The theory makes the 
following statements : 

1. For diffusion as a function of hydrostatic pressure, 
the diffusion coefficient is an exponential function of the 
volume strain. 

2. The rate of change of the diffusion coefficient with 
strain is related to the interatomic forces. The relation 
is explicit enough that the variation of the diffusion 
coefficient with pressure can be interpreted in terms of 
the interatomic potential-energy functions of the 
material. 

3. For diffusion under hydrostatic pressure, the 
activation volume can be calculated from the com- 
pressibility and the rate of change of the diffusion 
coefficient with volume strain. 

in every case for which data are available, these 
conclusions are in agreement with experiment. 

The general framework of the theory provides a 
basis for understanding the effect of strain on diffusion 
in terms of the atomic properties of the system and 
should provide a valuable tool for comparing diffusion 
rates for different states of strain, as well as for in- 
vestigating the mechanism of diffusion. 

APPENDIX 

The Vacancy Concentration Formula 

Consider a canonical ensemble containing X member 
systems, each system being a crystal containing N 
atoms and l vacancies. Let Ef be the jth energy level 
of a system containing l vacancies and let 0/ be the 
corresponding degeneracy. Then the number of systems 
containing l vacancies is 

Styexp (-Ef/kT) 

N,=X— , (Al) 

Z 

where Z is the total partition function for the ensemble. 

The number of vacancies in the ensemble is 

N, = E IN,, (A2) 

i 

and the number of atoms in the ensemble is 


Table II. Activation volumes for self-diffusion of various solids 
at 1 atmosphere calculated from Eq. (62). 


Solid 

Temperature, 

T 

<°K) 

Activation 

volume, 

AV\ 

(cc/g-atom) 

Sodium 

363 

12.3 

Phosphorous (white) 

314 

71.7 

Zinc (X) 

580 

3.0 

Zinc (||) 

580 

8.3 

Mercury (liquid) 

303 

0.62 

Gallium (liquid) 

303 

0.62 

Silver in silver chloride 

573 

10.3 

Silver in silver bromide 

573 

13.7 

Lead 

526 

13.9 

Lead 

574 

12.0 


N a =NZN,. (A3) 

i 

The atomic fraction of vacancies is given by n v =N,/ 
{N a +N v ). Since N v <&N a , n v is given by the ratio of 
Eqs. (A2) to (A3) to an excellent approximation, and 
therefore ' 

n.-ZWlfZQi), (A4) 

l I 

where Q, is defined by 

Q,= Z^i l eM-Ei l /kT). (A5) 

j 

Qi is the partition function of a system containing l 
vacancies. 
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Carrying out the division in (A4) and retaining only 
the leading term give 

«.= (l/WW^o, (A6) 

which is an excellent approximation, since the energy 
of formation of a vacancy is of the order of 1 ev, and 
therefore the higher terms in the series are very small. 
Qo is the partition function of a perfect crystal and Q i 
is the partition function of a crystal containing a 
vacancy. 

In the semiclassical approximation, 

~= (N+l) f f dpjdq, 

Q o J J j 

J...J g-MkT n dpjdqj, (A7) 

where p v =pv(p,q) is the energy of a crystal containing 



a vacancy and \po=Po (P a) is the energy of a perfect 
crystal. The integrations are carried out over all values 
of the momenta and coordinates pj and q,-. The factor 
(A+l) arises from the fact that N indistinguishable 
atoms can be placed in (A+l) numbered lattice sites 
in (iV+1) ways so that Qi is proportional to (A+l). 

Combining Eqs. (A6) and (A7) gives 


-/-/ e -M Pi.vmT JJ dpjdqj j 

J . . . J e-Mvi.^/kT j} dpjdqh ( A8 ) 


where unity has been neglected relative to N. 
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